The Fate of the Two-Magnon Bound State in the Heisenberg-Ising Antiferromagnet 
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The energy spectrum of the two-magnon bound states in the Heisenberg-Ising antiferromagnet 
on the square lattice are calculated using series expansion methods. The results confirm an earlier 
spin- wave prediction of Oguchi and Ishikawa, that the bound states vanish into the continuum before 
the isotropic Heisenberg limit is reached. 
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I. INTRODUCTION 

We consider the anisotropic antiferromagnetic Heisen- 
berg model on a bipartite lattice 



H 



<ij> 



<ij> 

Ho + xV 



(1) 



with S=l/2 spins interacting with their nearest neigh- 
bours. 

In the Ising hmit x = 0, the system is Neel ordered, 
with spins up {S^ ~ +1/2) on the A sublattice, let us 
say, and down {S^ = —1/2)) on the B sublattice. The 
1-particle 'magnon' excitations correspond to a single 
flipped spin on either the A sublattice, with total spin 
5^ = -1, or on the B sublattice (S*^ = +1). They have 
an excitation energy z J/2, where z is the lattice coordi- 
nation number. 

Separated two-particle excitations then have energy zJ 
in this limit; but a two-particle excitation on neighbour- 
ing A and B sites has energy only {z — 1) J, with total 
spin = 0, forming a 2-particle bound state. 

When one takes the isotropic limit a: —> 1, as is well 
known, the energy gap for the single magnon states van- 
ishes. The question then is, do the 2-particle bound 
states survive in this limit? 

It is well-known that two-magnon bound states exist 
for the isotropic Heisebnberg ferromapnet - see for exam- 
ple the textbook discussion by Mattis In recent years, 
there have also been considerable discussions of bound 
states in antiferromagnetic systems with frustration or 
anisotropy [1, ij. The question still remains, however, 
whether there might also be bound states in the simple, 
isotropic antiferromagnet. According to Mattis, this re- 
mained an "unanswered question" in 1965 (Ref. [l[, p. 
166). 

The question was investigated using spin-wave theory 
by Oguchi and Ishikawa [j| in 1973. They used lin- 
ear spin-wave theory with some fourth-order interaction 
terms included. They found that the two-magnon bound 
states merge into the continuum as one goes from the 
Ising limit to the isotropic limit, so that none survive. 



(a) (b) 



FIG. 1: Spin configurations on the linear chain in the Ising 
limit: a) Neel ground state; b) a single spinon state; c) two ad- 
jacent spinous, or a single magnon; c) two spinous separated 
by one site, or two adjacent magnons. 



This calculation neglects many higher-order effects which 
might be important, however, so it cannot be taken as 
definitive. Some numerical calculations involving multi- 
magnon states have also been carried out for the isotropic 
antiferromagnet on the square lattice. These include a 
studies of the spectral weights using series expansions Q 
and quantum Monte Carlo simulations @ , and a contin- 
uous unitary transformation (CUTS) study of the spec- 
trum Q. None of these works, however, have addressed 
the question of the bound states. 

In this paper we explore the fate of the two- 
magnon bound states for two particular cases, the one- 
dimensional chain and the two-dimensional square lat- 
tice. In one dimension, the model is exactly solvable, 
and the answer is already known: we simply review the 
results. In two dimensions, we use linked cluster meth- 
ods [1] to obtain series expansions in x up to order O(x^) 
for the bound-state energies, and extrapolate the results 
to a; = 1 using standard methods. In summary, our re- 
sults agree quite well with the spin- wave predictions 
The bound states disappear into the 2-particle continuum 
shortly before the isotropic limit x ^ 1 is reached. 



II. THE ONE-DIMENSIONAL CASE 

The one-dimensional model, i.e. the linear chain, is 
a special case. The model can be solved exactly us- 
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ing the Bethe ansatz [9|, and exact expressions for the 
low-lying spectrum have been obtained p^ . The in- 
dependent quasiparticles in this case are not 5^ = ±1 
magnons, but = ±1/2 'spinous', or domain walls. 
In the Ising limit x = 0, for instance, the ground state 
consists of alternating spins = ±1/2 on even or odd 
sites respectively, or the reverse (FigU^). A spinon or 
domain wall consists of a neighbouring pair of identi- 
cal spins in an otherwise alternating chain (FiglllD). For 
periodic or anti-periodic boundary conditions, spinons 
can only be created in pairs. Thus on an even lattice 
with periodic boundary conditions, say, the lowest-lying 
excitations above the ground state consist of a pair of 
spinons, which may be one lattice spacing apart (Figdt), 
two spacings apart (Fig[TJi), or any number apart. Fig. 
[TJ; then corresponds to a single spin-flip, and Fig. [TJi to 
two neighbouring spin-flips, and so on. The excitation 
energy in this limit, however, is AE' = J, independent of 
the spacing. 

Johnson, Krinsky and McCoy [lO| have obtained an 
exact expression for the low-lying excitation energy at 
all x: 

AE = ^^Lfi^ [(1 _ cOS^ + (1 - fc2 ^Qg2 ^^)l/2] 

TT cosh l> 

(2) 

where cosh;^ — 1/x, K is a. complete elliptic integral of 
the first kind whose nome is q = E~'^ and modulus fc^, 
and < qi,q2 < TT are two free parameters corresponding 
to the momenta of the two domain walls in each sector. 
This spectrum is the same in the S"^ = ±1 and the (dou- 
bly degenerate) = sectors, so that the two S — 1/2 
spinons combine to form a degenerate S' = 1 triplet and 
S = singlet, at all couplings x in the range. So in 
this case, the state with two neighbouring flipped spins 
merely forms part of the 2-spinon continuum. 

III. THE SQUARE LATTICE 

As the prime example of a bipartite system in two di- 
mensions, we consider the model on a square lattice. In 



this case no exact solution is known, and so we have 
calculated numerical estimates for the energy of the 2- 
magnon state using series methods 0. We perform an 
Ising expansion, taking the Ising Hamiltonian Hq in equa- 
tion (P) as our unperturbed starting point, when the 
bound state consists simply of a pair of flipped spins on 
neighbouring sites, as discussed above. A perturbation 
series expansion in x is then calculated for the bound 
state energy, with V in equation ([J) as the perturbation 
operator. As a technical point, we note that the bound 
state lies in the same sector as the ground state, and 
hence a 'multiblock' diagonalization algorithm [s^ must 
be employed. 

Note also that there is one bound state configuration 
for each lattice bond, making a total of four times as 
many configurations as for either of the single-magnon 
states. Correspondingly, we obtain results for four dif- 




FIG. 2: Cuts through the Brillouin zone Fi — r4 correspond- 
ing to the four two-magnon bound states. 



ferent paths Fi, • • ■ r4 in the Brillouin zone, as shown in 
Fig. [21 whereas only the Fi mode is independent for the 
single-magnon state. 

The calculations have been carried out through 0{x^). 
Since only even-order terms appear, this corresponds to 
only five series coefficients at any fixed momentum. The 
leading order terms in the dispersion relation for the 
bound state excitation energy are: 



^ + ^ cos Y cos Y + ^(cos + cos ky) 



-— cos kx cos ky H (cos 2fc„ + cos 2ky] 

6 24 



1 , k^ 3 
-(cos — cos -ky 
6^ 2 2 ^ 



cos -kx cos — I 
2 2 ' 



(3) 



The complete series coefficients at selected momenta are 
listed in Table H 

Estimates of the bound-state energy as a function of x 
were now obtained using Fade approximants to extrapo- 
late the series. Since the number of coefficients is small. 



the accuracy of the extrapolation is also low. At smaller 
values of x, nevertheless, quite good estimates are possi- 
ble. For example. Figure O shows dispersion relations for 
the four bound-stste modes at a:; = 0.8, as compared with 
the lower bound of the 2-particle continuum. It can be 
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seen that all four modes remain bound below the inter- 
mediate plateau of the continuum, and only near k = 
do three out of the four modes merge into the continuum. 
The first mode appears to remain bound at all momenta. 

At X — 1.0 it is a different story, as shown in Figure 
m In short, it appears that none of the four modes re- 
main bound at any momentum. The nominal error bars 
are much larger in this case, firstly because x is larger, 
but also because we may expect some sort of singular be- 
haviour where the bound state merges with the contin- 
uum. The clearest picture is obtained for the first mode, 
where near k = the estimates lie on the lower edge of 
the continuum, within errors. There is an apparent lev- 
elling off at very small momenta, but this is a common 
feature of Fade extrapolations in Heisenberg-type mod- 
els; if a Huse transform ^llj was performed before the 
extrapolation, something much closer to the continuum 
behaviour would be expected. At larger momenta, the 
estimates generally lie well above the continuum lower 
bound. 

These results may be compared with the spin-wave 
results of Oguchi and Ishikawa who calculated the 
bound-state energies as functions of x for two specific mo- 
menta. At k = (0,0), they found that the bound states 
merged with the continuum in the range 0.7 < x < 0.95; 
we find only one bound state remaining at a: = 0.8, in 
agreement with their results. At k = (7r/2,7r/2), they 
found the merger to occur at somewhat larger values 
0.96 < X < 0.98; we find that all four states remain 
bound at X = 0.8, but have vanished by a; = 1.0, once 
more in agreement with their results. Of course, the ac- 



tual values for the energies have changed as higher-order 
terms are added in, but not by a large amount. 



IV. CONCLUSIONS 

We have used series expansion methods to calculate the 
energy of two-magnon bound states in the anisotropic 
Heisenberg-Ising antiferromagnet on the square lattice. 
We find that the bound states do not survive in the 
isotropic limit, in agreement with the spin- wave predic- 
tions of Oguchi and Ishikawa [1]. There are no bound 
states in the linear chain model either. Hence one may 
extrapolate that there will be no bound states in the 
isotropic Heisenberg antiferromagnet on any bipartite 
lattice, in sharp distinction to the ferromagnetic case. 
Oguchi and Ishikawa [4] have given some qualitative ar- 
guments why this might be so. 
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FIG. 3: Spectrum of the two-particle bound states Fi — r4 at a; = 0.8. The dashed Unes mark the lower edge of the 2-particle 
continuum. 



TABLE I: Ising expansion series coefficients in powers of x 
for the excitation energy A£J of the 2-particle bound state at 
selected momenta k = {kx,ky). 



[kx, ky) 


(0,0) 


(V2,0) 


(^,0) 


(7r/2,7r/2) 


(7r,7r) 





3.00000000000000 


3.00000000000000 


3.00000000000000 


3.00000000000000 


3.00000000000000 


2 


-1.66666666666667 


-0.804737854124365 


-0. 166666666666666 


-0.333333333333334 


-0.333333333333333 


4 


0.299074074074084 


0.412245678974668 


-0.0712962962962898 


0.403819444444453 


-0.727546296296289 


6 


-2.21500154321004 


-0.955337158046685 


-0.371659167631400 


-0.517003970550667 


-0.576196887860336 


8 


5.95191488216504 


2.06716252423366 


0.237274308271624 


0.738026171850790 


-0.604532331023972 
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FIG. 4: Spectrum of the two-particle bound states Fi — r4 at a; = 1.0. The dashed Unes mark the lower edge of the 2-particle 
continuum. 



